Analog of Astrophysical Magnetorotational Instability in a Couette- Taylor Flow of 

Polymer Fluids 
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We report experimental observation of an instability in a Couette- Taylor flow of a polymer fluid in 
a thin gap between two coaxially rotating cylinders in a regime where their angular velocity decreases 
with the radius while the specific angular momentum increases with the radius. In the considered 
regime, neither the inertial Rayleigh instability nor the purely elastic instability are possible. We 
propose that the observed "elasto-rotational" instability is an analog of the magnetorotational in- 
stability which plays a fundamental role in astrophysical Keplerian accretion disks. 
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Introduction. — Accretion is a fundamental process in 
astrophysics by which protostellar objects and stars are 
formed. Due to gravity, the interstellar gas collapses into 
thin disks differentially rotating around accreting bodies. 
In order for the gas to further fall onto the central object, 
angular momentum has to be transported out of the sys- 
tem. Since the molecular viscosity in disks is very small, 
the laminar Keplerian disk cannot loose its angular mo- 
mentum on astrophysically reasonable timescales. The 
need for much larger, possibly turbulent angular momen- 
tum transport was identified in Lynden-Bell [l[ , Shakura 
& Syunyaev 0] , although it remained unclear what could 
make a hydrodynamically stable Keplerian flow turbu- 
lent. In the 1990's it was realized that a weak magnetic 
field existing in accretion disks leads to a quickly grow- 
ing instability rendering the disks turbulent. This is the 
magneto- rotational instability (MRI) originally derived 
by Velikhov |i| and Chandrasekhar jj] and rediscovered 
in astrophysical context by Balbus & Hawley [j| . 

During the last decade, considerable progress has been 
made in understanding the effects of such instability on 
differentially rotating flows. A sizable amount of analytic 
work was devoted to linear analysis of instability thresh- 
olds [e.g., 0, 0, HI- Extensive numerical simulations of 
the nonlinear stage of the magnetorotational instability 
have also been performed, e.g., [I, E3, U, 12, EH, how- 
ever, it still remains a challenge to address the ranges 
of scales relevant for real astrophysical systems. Labo- 
ratory investigations of the magneto-rotational instabil- 
ity in li quid -metal experiments have been proposed [e.g., 
El El, til E3, EH and conducted 0, Ej. However, 
large resistivity of liquid metals complicates unambigu- 
ous laboratory study of the magnetorotational instabil- 
ity. Recently, however, the instability was observed in a 
Couette- Taylor liquid metal experiment where the helical 
rather than "standard" axial magnetic field was applied 
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by external coils [21(. The relevance of this setting for 
Keplerian accretion disks is discussed in [22| ■ 

In the present paper we report a laboratory observa- 
tion of an analog of astrophysical magnetorotational in- 
stability in an experiment using visco-elastic solutions of 
high molecular weight polymers. In a certain range of pa- 
rameters, the dynamic equations describing visco-elastic 
polymer fluids arc identical to the magnetohydrodynamic 
equations describing conducting fluids or plasmas. This 
opens a way to investigate the fundamental astrophysical 
instability in a simple laboratory setting. 

To explain the physics of the instability, consider two 
fluid elements rotating at different orbits and connected 
by an elastic string (a magnetic field in accretion disks, a 
polymer in our experiment). The inner element rotates 
faster, therefore it is pulled back by the string. As a 
result, it loses its angular momentum and falls closer to 
the center. The outer fluid element is pulled forward, 
gains angular momentum, and goes to a larger orbit. The 
fluid elements thus move apart stretching the string even 
more, leading to the instability. 

Our interest to the problem was motivated by analytic 
work of Ogilvie & Proctor [23| elucidating the analogy 
between instabilities in Couette- Taylor flows of magnetic 
and polymer fluids. [When our experiment was in 
progress, we become aware of the new paper by Ogilvie 
& Potter 24( where this analogy is developed in more 
detail.] The experiment provides an intriguing possi- 
bility to investigate the regime of "elasto-rotational" 
instability and resulting "elasto-rotational" turbulence 
in non-Newtonian polymer fluids; such regimes have not 
been experimentally studied before. 

Magnetohydrodynamics and polymer fluid dynamics. — 
The dynamics of a conducting fluid is described by the 
set of magnetohydrodynamic (MHD) equations: 

<9 f v + (v • V)v = -Vp + (B • V)B + ??Av + F, (1) 
d t B + (v • V)B - (B • V)v = 77 M AB, (2) 

where v(x, t) is the velocity field, B(x, t) is the magnetic 
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field normalized by ^/AiFp, p is pressure which includes the 
magnetic part, r\ is fluid viscosity, and r\M resistivity. We 
assume that the fluid is incompressible and the density is 
constant, say p = 1. The external force F(x, t) represents 
mechanisms driving the flow. We also assume cylindri- 
cal geometry with the coordinates r, 9, and z, where 
the steady state is described by the azimuthal velocity 
field vg(r). For the gravitational force, F r oc — 1/r 2 , the 
velocity field has the Keplerian profile vg(r) oc r -1 / 2 . 

As follows from Eq. (JTJ , the back reaction of the mag- 
netic field on the flow is described by the Maxwell elec- 
tromagnetic stress tensor Tli = B % BK The evolution 
equation for this tensor is derived from Eq. Ip]). where 
we neglect small resistivity tjm- 



M + (v • V)T« - T% - I* Vl^ = 0. 



(3) 



The tensor TVL = B l B^ obeying this equation is said to 
be "frozen" into the flow. In the case of polymer fluids, 
the polymer stress tensor Tp frozen into the flow should 
obey the same equation. 

In contrast with a magnetic fluid, there is no exact 
equation describing the dynamics of polymer solutions. 
However, in the case when a solution is dilute, one can 
formulate the constitutive equations based on general 
principles of fluid dynamics [3, [2(| . The stress tensor 
can be represented as a linear sum of the viscous stress 
of the solvent T ij = ^{V^ + V 3 v l ) and the stress T l J 
contributed by the polymer. The polymer contribution 
should generally obey the equation: 



1 p 



(4) 



Here r is the relaxation time of the fluid element, which 
is related to polymer elasticity, and D t denotes the con- 
vective derivative, as in Eq. (|3"|), 



D t rj = dtTp 1 + (v • v)rj 



TUV,v j . 



(5) 



If the relaxation time r is very large compared to a 
characteristic time of the flow, the second term in the 
left hand side of Eq. {4| dominates, and the stress is 
advected by the fluid. In the other limit, r — > 0, the 
polymer is not frozen into the fluid - it rapidly relaxes 
to its non-stretched equilibrium configuration and con- 
tributes to fluid viscosity, Tp = ^(V^w- 7 + Vji/). The 
dynamics of the velocity field in Eq. ([5]) is given by the 
standard Navier-Stokes equation of motion: 



work, we change the variable: Tp — > Tp = Tp + ^-S 1 ^ , 
where ft 13 is the Kronecker delta. The momentum equa- 
tions flTJ, © for magnetic and polymer fluids now have 
identical forms: 

<9 4 v + (v ■ V)v = -Vp M + V ■ T M + ijAv + F, (7) 
a t v+(v V)v = -Vpp + V-Tp+r/Av + F, (8) 

where the pressure terms ensure incompressibility of the 
flows. The dynamic equations for the stress tensors are: 



^ + (v-V)^ 



■ M 



= r lM [B l V 2 B 1 +(V 2 B l )S J ], 



d t rj + (v • v)rj 



T^Viv 1 



1 r,H,4 

T 



(9) 



(10) 



These equations are identical except for the dissipa- 
tion terms - the magnetic field diffuses while the poly- 
mer stress relaxes. However, if the magnetic Reynolds 
number Rm ~ flcP/rjM and the Weissenberg number 
Wi = r|9fi/9(lnr)| ~ rfl are large (CI being the an- 
gular velocity and d the gap between cylinders), one can 
neglect the dissipation terms. 

Denote R± and i?2 as the inner and outer radii, respec- 
tively. When the gap is narrow, d/R <§C 1, the basic veloc- 
ity flow can be represented as v e (r) = v (Ri) + (r — Ri)Cl. 
The corresponding stationary solution of Eq. (|10[) in co- 
ordinates (r, 9, z) then has the form [231 ] : 



T 



1 -Wi 
-Wi 2Wi 2 + l 
1 



(11) 



There is no exact correspondence of the tensor (fTTj) 
to the magnetic tensor Tm, since (fTTj) cannot be repre- 
sented as a product of two vector fields. However, one 
can introduce a set of three auxiliary fields, B 1; B 2 and 
B 3 such that 
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B\B{ 



BlBl 



(12) 



In this representation, Bi and B 2 have radial and az- 
imuthal components, while B3 is purely axial |23j . 



Bi, 



(?) 



1/2 



Wi ± (Wi 




1) 



1/2 



(13) 



d t v + (v ■ V)v = -Vp + V • T P + 77AV - 



(6) 



The system (HJ), ([5]) and ^ presents the so-called B 
model of Oldroyd [26|, a constitutive system for dilute 
polymer solutions. 

The Ogilvie- Proctor model. — The analogy of MHD and 
polymer fluid instabilities in the Couette- Taylor regime 
was elucidated by Ogilvie & Proctor [23| . Following their 



B, 



1/2 



(14) 



A general analysis of the instability requires expan- 
sion of the nonlinear equations JHJ) and (I10| in small 
deviations from the basic flow. Depending on what de- 
viations are considered, different "magnetic fields" play 
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dominant roles. If one assumes that the perturbations 
are axisymmetric (kg — 0), and their wavevectors obey 
k z ^S> k r , then the azimuthal and radial fields B 12 are 
not relevant for the instability, and the dominant role is 
played by the axial field B3, in direct analogy with the 
corresponding magnetorotational instability. 

Experiment. — In our experiment, a polymer fluid fills 
the gap between two coaxial cylinders rotating in the 
same direction with different angular velocities. The gap 
is narrow, and the cylinders are driven by the same mo- 
tor with two different gears to approximate the Keplerian 
velocity profile O(r) oc r -3 / 2 . (In fact, any profile where 
the outer cylinder rotates slower than the inner one, but 
the specific angular momentum of the outer cylinder is 
larger than that of the inner one is suitable for the con- 
sidered instability.) 

The outer diameter of the inner cylinder is 14" , the in- 
ner diameter of the outer cylinder is 15.5" , and the height 
of the cylinders is 2'. The outer cylinder is transpar- 
ent and the flow is visualized by adding a small amount 
of highly reflecting Kalliroscope particles. The angular 
velocity of rotation can reach 40rad/s, which, for the 
Keplerian velocity profile translates into a shearing rate 
7 = dvg/dr — vg/r « 60s -1 . 

For the polymer fluid we choose an aqueous solution 
of high molecular weight Polyethylene Oxide (Mw ~ 
7,000,000) obtained from DOW Chemical. The exper- 
iments were conducted at ambient temperature of 20 C, 
although the temperature was not precisely controlled. 
First, we checked that in the studied range of angular 
velocities, the hydrodynamic flow without polymer addi- 
tives was stable. We then performed a series of experi- 
ments with different concentrations of PolyOx. In each 
experiment we gradually increased the rotation velocity 
to obtain the instability threshold. 

No instability was observed for concentrations less 
than about 0.2% by weight, rather, at very high rota- 
tion rates turbulence set up at the ends of the cylinders 
where the Keplerian profile is broken, and propagated 
over the whole cylinder. At higher concentrations, how- 
ever, the instability did appear. At the polymer concen- 
tration of 0.25 %, the most unstable mode was a spiral 
v(k) oc exp(ik z z + imO) with azimuthal wave number 
m = ±1 and the axial wavelength A 2 « 20 mm. Due to 
symmetry, the spirals winding up and down are equally 
probable. In different runs, the flow therefore sponta- 
neously broke into regions of m = 1 and m — — 1, as 
e.g. in Fig (p}, left panel. The threshold for this instabil- 
ity was about 7 C (0.25) ~ 5.6 s -1 . As the concentration 
was increased further, the most unstable mode became 
axisymmetric with the wavelength X z rs 30 mm. The 
threshold for this instability was y c (0.5) » 7.3 s -1 . The 
result for 0.5 % solution is shown in Fig ([I]), right panel. 
In both cases the instability was detectable by eye and 
the pattern was captured with a generic digital camera. 




FIG. 1: Elasto-rotational instability in visco-elastic Couette- 
Taylor flow. Left: Most unstable mode is m = ±1. Right: 
Most unstable mode is axisymmetric. 

To argue that the observed instability is analogous to 
the magnetorotational instability, we performed another 
series of experiments. This time the gears at the cylinders 
were chosen to set up either quasi-Keplerian O(r) oc r -1 ' 3 
or "anti-Keplerian" f2(r) oc r 13 profiles. This is done 
to exclude the so-called elastic instability that can exist 
even for small Reynolds numbers (Re = fid 2 jr\p) as long 
as the Weissenberg number exceeds a certain threshold 
[27I] , The elastic instability takes energy from the elastic 
energy of the flow, and should not essentially depend on 
the sign of dfl/dr, while for the MRI the sign of dfl/dr 
is crucial. 

In the "anti-Keplerian" case, the only instability 
that could exist is purely elastic instability. For both 
considered concentrations of PolyOx, we observed the 
instability in the quasi-Keplerian case (analogous to the 
instability in the Keplerian case), however, we did not 
observed any instability in the anti-Keplerian case, even 
when we increased the shearing rates to three times as 
high as in the Keplerian counterpart. This indicates 
that in our Keplerian case the observed instability is 
driven by inertia and takes its energy from the kinetic 
energy of the flow. We therefore propose that the 
observed "elasto-rotational" instability is analogous to 
the magnetorotational instability. 

Discussion. — Certain support for our observations is 
provided by the Ogilvie-Proctor consideration outlined 
in previous sections. We should be cautioned, however, 
that this model is to some extent phenomenological. It is 
known that the polymers fluid viscosity r\ v and relaxation 
time t are not constants, but they strongly decrease as 
the shearing rate increases beyond 7T ~ I (the effect of 
shear-thinning). Besides, a model with single relaxation 
time is often not adequate, and one needs to introduce a 
series of relaxation times describing the relation between 
the shear rate and stress tensor. More essential, however, 
is the fact that in our case the polymer solution cannot 
be considered dilute for the used polymer concentrations. 
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The application of the theory is therefore limited. 

We however found that the theory is in reasonable 
agreement with the experiment if one substitutes the ex- 
perimentally measured value for unknown viscosity r] p . 
Indeed, consider the case of the axisymmetric instability 
observed at polymer concentration 0.5%. We measured 
the shear viscosity of the flow rj p at the obtained critical 
shearing rate 7 C by measuring the viscosity in a Bohlin 
rheometer using a system of two coaxial cylinders - a 
scaled down copy of the experimental set-up. The inner 
cylinder was stationary, and the outer cylinder rotated 
steadily so that the shear rate in the gap matched the 
shear rate in the experiment. The measured viscosity 
was ?7 P (0.5) ~ 187mPa ■ s. As for the relaxation time, it 
can be found for very low shearing rates using oscillation 
measurements, giving the value of order to(0.5) ~ 1.4s. 
The relaxation time is however strongly shear-thinned at 
the experimental shear rate, so its precise value is difficult 
to evaluate. 

We now substitute the experimental values of the criti- 
cal shearing rate j c , viscosity r/ p , and the wave vector k z , 
into the linearized Oldroid-B equations. Such linearized 
equations are derived in the limit of small but nonva- 
nishing d/R in (27l |: they are bulky and not presented 
here. We solved these equations numerically. The solu- 
tion confirms that the axisymmetric instability with the 
observed parameters indeed exists if the fluid relaxation 
time is r ~ 0.6s. This is a reasonable number if shear- 
thinning is taken into account. With this relaxation time 
we estimate that the instability occurs at Wi = -yr ~ 4.4 
and Re — 7c? 2 /rjp ~ 14. Moreover, when we numerically 
switched to the anti-Keplerian profile by inverting the 
sign of the shearing rate 7, the instability disappeared, 
which agrees with the experiment. 

We note a useful fact that in the axisymmetric case 
(m = 0), the instability threshold involves only B3 = 
\A?p/ r - I n the kinetic theory of dilute polymer solu- 
tions the ratio rjp/r is constant and proportional only to 
polymer concentration [e.g.,[28|. Therefore the "imposed 
magnetic field" B 3 is stable even when both polymer vis- 
cosity and relaxation time are shear-thinned by the flow. 
In the non- axisymmetric case (m ^ 0), the instability 
also depends on the azimuthal fields B1.2. In principle, 
it may be possible to design an experiment where such 
azimuthal field dominates, in even closer analogy with 
real accretion disks. 

In conclusion, based on our results we propose that 
the analog of magnetorotational instability can be exper- 
imentally studied in visco-elastic flows of polymer fluids. 
For a more quantitative analysis and for direct compar- 
ison with the theory, different polymer solutions whose 
viscosities are not strongly shear-thinned should be used, 
the so-called Boger fluids [2§|. This work is in progress. 
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